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Abstract 

The recent Tevatron result on the top quark forward-backward asymmetry, which deviates from 
its standard model prediction by 3.4cr, has prompted many authors to build new models to account 
for this anomaly. Among the various proposals, we find that those mechanisms which produce ti 
via t- or ii-channel can have a strong correlation to the rare B decays. We demonstrate this link 
by studying a model with a new charged gauge boson, W' . In terms of the current measurements 
on B — > nK decays, we conclude that the branching ratio for B~ — > ir~ K° is affected most by 
the new effects. Furthermore, using the world average branching ratio for the exclusive B decays 
at 2cj level, we discuss the allowed values for the new parameters. Finally, we point out that the 
influence of the new physics effects on the direct CP asymmetry in B decays is insignificant. 
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Recently, the forward-backward asymmetry (FBA) in top-pair production is measured 
by the D0 |l| and CDF j^j Collaborations in pp collisions at y/s = 1.96 TeV. With an 
integrated luminosity of 5.3 fb _1 , CDF further finds that the large top quark FBA occurs 
at large tt rapidity difference (Ay) and invariant mass of tt (M«) [3J|: 

A a (\y\ < 1.0) = 0.026 ±0.118 [0.039 ±0.006], 
A*(\y\ > 1.0) = 0.611 ±0.256 [0.123 ±0.008], 



A a (M tI < 450 GeV) = -0.116 ± 0.153 [0.040 ± 0.006] 
A tl {M t i > 450 GeV) = 0.475 ±0.114 [0.088 ± 0.013] , 



(1) 



where the value in the square brackets denotes the theoretical result calculated by Monte 
Carlo program (MCFM) |4| at next-to- leading order (NLO). The inconsistency between data 
and standard model (SM) predictions displayed in Eq. (CQ) indicates that new physics effects 
may be at play and hence leading to this anomalous result. 

In the wake of this 3.4<r deviation of the observed value of the top quark FBA from the 
SM predicted one, several possible solutions have been proposed and studied by authors 



in Refs. |s-l30l. I32l44lj . Among these new mechanisms, we see that a potentially interesting 
correlation to B meson physics may arise. In particular, the new interaction that contributes 
to ti production by t- or u-channel will also contribute to rare B decays, e.g. B — > 7i(ir, K). 
This is achieved thr oug h box diagrams like those sketched in Fig. [TJ where X could be a 



a tnrougn dox diagrams like tnose sketcnea m t lg. LU wnere i\ 
colored vector fl, H, 0], Z' |§j, W 0, 30, 32] or colored scalar bosons [ll, [l2|, llJ, |35|, 
while q' = d or u depending on the charge of the X particle. It is known that with the 
enormous B meson production at LHCb, BaBar, Belle and Tevatron, the errors in the 
measured branching ratios (BRs) for B — > nK decays have now reached percent level. As 
a result of this high accuracy, it is interesting to investigate how strong is the correlation 
between rare B decays and top quark FBA. In addition, one can study the influence on the 
associated CP asymmetry (CPA) in such new physics setups. We present our analyzes on 
both issues in the following. 

In order to illustrate the impact of the new physics (which leads to a large top-pair FBA 
as measured) on the low energy B sector, in this paper, we shall focus on the case with 
a new W interaction, namely, the t-d-W coupling. A similar discussion could be applied 
to other new interactions in the Z' colored vector or colored scalar models. Since we are 



studying the impact of new physics on B decays, the detailed analysis of top-pair production 



FIG. 1: Feynman diagrams for b — > q'q'q with q' = d or u and q = d,s, where X stands for a 
generic new particle. 



could be referred to those papers listed in the references. We start by writing the relevant 



interaction as 



30] 



£ = -i^t^WlP R d + h.c , (2) 

where g' 2 is the new gauge coupling which is regarded as a free parameter. By including 
the contributions of W and the associated Goldstone bosons, the Hamiltonian for b — > qdd 
(q = d, s) shown in Fig. [I] is given by 

Uw = ~^ V ^k{^^) I(xw,x t )d(l- l6 )bq(l + ls )d, 

= -9z Vtb V t * q C w/ d(l - lb )bq{\ + l5 )d (3) 

where Cw> stands for the new Wilson coefficient at electroweak scale, xw = m^/m^,, 
x t = mt/rriw' and I(x w ,Xt) = (1 + 2(x w + Xt)h{xwi%t) with 

-1 rzi 



h(a, b) = / dz\ I dz 2 
Jo Jo 



~2 



z 2 



(1 — a)zx — (a — b)z 2 ' 

h(a,b) = I dz x I dz 2 - ^ 2 — — ^ . (4) 

Jo Jo (1 - (1 - a)Zi -{a- b)z 2 y 

We note that there are only two new free parameters in Cy/ 1 - To illustrate the dependence 
on g' 2 and the mass of W, we plot the contours for Cy/ 1 as a function of g' 2 and mw' 111 
Fig. [2J The number on each curve in Fig. [2] denotes the value of Cw> in units of 10 -2 . For 
instance, with g' 2 = 3 and m W i = 700 GeV, we get Cw 1 = 0.85 x 10 -2 . 

Notice that if we perform Fierz transformation on the fields in Eq. ([3]), the four- Fermi 
operator could be transformed into d a 7^(1+75)^/3^7^(1— 75)^, where a, j3 denote the color 



indices. This operator is the same as the Oq operator in the SM 



43| and its corresponding 



Wilson coefficient at mw scale is Cq ~ —0.2 x 10~ 2 . According to our brief analysis, we see 
that the new interaction has the same behavior as Oq, and so, its contribution should be 
similar to that dictated by C 6 . 
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FIG. 2: Contours for Cy/i as a function of g' 2 and mw, where the number on each curve denotes 
the value of Cw in units of 10~ 2 . 

Since the new four- Fermi interactions in Eq. fl3]) are induced by loop diagrams, we expect 
that the influence of the new effects on the tree dominant B — > nn decays should be small. 
Therefore, we only concentrate on the B — > 7rK decays. This is because at quark level, 
they are associated with the b — » sdd process which is dominated by gluonic penguins in 
the SM. In order to comprehend how the new effects contribute to each decay mode in the 
B — > txK processes, we draw the flavor diagrams in Fig. [3j where figures (a) and (b) denote 
the color-allowed and -suppressed emission diagrams respectively, while figure (c) represents 
the annihilation diagram. Importantly, we note that the decay B~~ — > tt~K°, B~ — > tt°K~ 
and B d — > n + K~ are influenced by Fig. E|a), (b) and (c) respectively, whereas B d — > 7i°K° 
involves all three figures. 

It is well-known that the calculations of a nonleptonic exclusive decay include factorizable 
and nonfactorizable parts. Owing to the color flow, the latter is always color-suppressed and 
smaller than the former. Hence, to simplify the analysis, we only calculate the factorizable 
contributions from the W exchange when displaying the influence of the new physics effects. 
However, we keep both parts when accounting for the SM sector. In order to study the 
nonleptonic exclusive B decays, we parametrize the relevant decay constants, transition and 
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FIG. 3: Flavor diagrams for B — > ttK decays: (a) color-allowed emission, (b) color-suppressed 
emission, and (c) annihilation diagrams. 

time-like form factors as 



(0|g 2 7 M 759i|M> = if M P^ 
(M|g l7 M 75g 2 |0) = -if M P^ 



(M'\q 2 <fq 1 \M) = /f M \q 2 ) ( P» - 1_V ) + {q 



P-q 



{M.M^qM = g^(Q 2 ) ( ? - ) +g^(Q 



q 

Q-q 



vMM' (J2- 



P ■ q 



MiM 2 / n 2\ Q • q 

Q2 J > *0 W ) Q2 



Q 11 



(5) 



with P = P M + Pm', q = Pm ~ Pm>, Q = Phh + Pm 2 and q = P Ml - P Ah . Applying the 
equations of motion, we then get 



m 



M 



(0|g 275 gi|M) = (M|g l75 g 2 |0) = -if M - 



(M'\q 2qi \M) 
{M 1 M 2 \q 1 q 2 \0) 



M 2 - M' 2 

M 2 - M 2 

m m - m, 



/o MM '(g 2 ), 
9 A o hAh (Q 2 ) 



(6) 



"qi ""<}2 

Using the above form factors, the decay amplitude with W exchange is formulated by 



_ G F V tb V t * 



A w '(nK) 



V2 



(7) 



for B — > Ti K°, B>d — > tt + K and B — > 7r°K processes, while 



A W '(7T°K ) 



V2 



A (ir~K ) + A (ir K~ 



V2 



A w '(tt + K- 



(8) 
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for Bd — y 7T°K , where 

2 2 2 
iir-K" - Ik : Jo \ m K) > 

m s + m d m b - m u 

m| m\-ml Kn 2 

Ctt+ZT- - JB j P \ m B) ) 

m b + m u m s - m d 

2 2 

C.ok- - A 2V5jVe Jo ("O • (9) 

Clearly, once the first three decays have been determined, the last one will be fixed. Since 
m K it ^ m Bi it i s a good approximation to take fo™ (m'jcM) ~ fa* (0)- Moreover, the 
transition form factor /(f M (0) an d time-like form factor g^ lM2 (rn B ) are associated with SM 



QCD, and for calculating them, we employ the perturbative QCD approach (PQCD) 44]. 
The resulting formulae are summarized in Appendix |A] We note that in Eq.(J7J), Cw'(^b) is 
the new Wilson coefficient at /i& = 0(mj,) scale. 

Now that the necessary formulations have been built up, we list the SM theoretical inputs 
required for our subsequent analysis in Table [H where V ub = \V u b\e~ l< ^ 3 . Using these values 
and formulations in the Appendix, the transition and time-like form factor are found to be 

fB"(0) « 0.24, 
f BK (0) « 0.35, 
gl K {ml) « 0.15e iL1 . (10) 

For estimating CV'(/x&), we adopt the relation: Cw'(^b) ~ (®s(Hb) / ®s( m w)) d Cw'( m w) with 
d = N C (N 2 — l)/(lliV c — 2nf), N c = 3 and rif being the number of effective quark flavors 



491 ] . In order to combine the SM results with the effects of W exchange, we write the total 

TABLE I: The values of theoretical inputs. 



\V ub \ 


03 (deg) 


Vts 


m% [GeV] 


m° K [GeV] 


3.9 x 1(T 3 


65 


-0.041 


1.34 


1.70 


m d{u) [MeV] 


m s [MeV] 


U [MeV] 


fx [MeV] 


fs [MeV] 


4.5 


100 


130 


160 


190 



decay amplitude for B — y irK as 

A(ttK) = A sm (ttK) + A w '(ttK) . (11) 
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As usual, the SM contributions could be expressed by [46] 



A SM (ir~K ) 
A SM (n + K~) 



P. 



GpVtbVfg 

G F VtbV* 8 _ G F V ub V* s 



V2 
G F V tb V* 



V2 



T . 



V2 



{P + Pew)~ GfVu ^ :s {T + C) 



V2A SM (ir°K°) = g^S ( P _ Pj 



y/2 
GpVubV* 



V2 - * W 

In Eq. (112]) . P, Pew j T and C denote, respectively, the contributions from gluonic pen- 
guin, electroweak penguin, color-allowed tree and color-suppressed tree. To deal with these 
hadronic effects, we quote the results calculated by PQCD, where the values in the SM are 



C. 



(12) 



given by 



47] 



P = 0.15e 



-i0.24 



P EW = 0.018e 



i0.44 



1.05e 



iO.l 



C = 0.270e 



-il.3 



(13) 



Consequently, the resulting SM 
with the experimental data 48] 



predicted values of the BRs and CPAs are presented alon^ 
in Table M 



TABLE II: SM results [47[ and experimental data [481 ] for BRs and CPAs of B — > ttK decays. 



Decay 


Br -> vr-i^° 


5 d -> vr+K- 




5 d -> 


6 SM [10- 6 ] 


23.5 


20.46 


13.25 


9.16 




23.1 ± 1.0 


19.4 ±0.6 


12.9 ±0.6 


9.5 ±0.5 


A™[%] 





-10.18 


-0.95 


-6.29 


A E c x P p [%] 


0.9 ±2.5 


-9-811;? 


5 ±2.5 


-1 ± 10 



To display the influence of the W effects on BRs and CPAs in B decays, we consider the 
two quantities Rb and Rcp which represent, respectively, the ratios of BR and CPA in the 
W exchange to those in the SM: 

B W '(B -»• tiK) 



Rb 
Rcp 



B SM (B -> vr K) 
A^(B^nK) 



(14) 
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Because there is only two new free parameters in the W-mediated effects, we use contour 
plots as a function of g 2 and mw 1 to display the deviation from the SM predictions. With 
the SM inputs discussed earlier, the contours for i?^ as a function of g' 2 and m^i are shown 
in Fig. HI where (a)-(d) correspond to the decays B~ — y tt~K°, Bd — y ir + K~ , B~ — y ir°K~ 
and Ed — y 7r°K°, respectively. From this, we see that W- mediated effects have a significant 
contribution to the BRs of the tt~K°, tt°K~ and tt°K modes, whereas the effects are small 
for the n + K~ mode. The main reason for this can be traced to the fact that the associated 
topology for Bd — y tt + K~ decay is of annihilation type, whose corresponding hadronic effects 
are usually smaller than those arising from the emission diagrams. 

In addition, if we compare the three decay modes -k~K°, tt°K~ and -k°K°, one finds 
that the influence on ir°K~ mode is smaller than the other two because Fig. EJb) is color- 
suppressed. From our numerical results, we conclude that B~ — y -k~K° decay is the most 
sensitive to the W^'-mediated effects. In order to display the constraint from the measured 
data, we plot the allowed range for g' 2 and m^y in Fig. |5l where we have adopted the world 



average BR for B~ -> ti~K° with la errors, (2.31 ±0.20) x 10~ 5 [31|, and the PQCD results 
with the uncertainties are taken as \P\ = 0.15+° f 3 and arg(P) = -0.24^ Q. From these, 
we find that only those values of \P\ approaching to the central value calculated by PQCD 



contribute to the range 400 < mw' < 800 and the allowed m 



v 



i from other values of \P\ is 



30] . Based on the result, we 



below 200 GeV, which is disfavored by the analysis in Ref. 
observe that g' 2 > 2.5 is excluded when mw < 700 GeV. Such constraint on the parameter 
space highlights the tension between the need to choose a large enough coupling (g' 2 — 3 



when m' w ~ 550 GeV 32J) for this type of model to explain the top anomaly and satisfying 
the limits implied by rare B processes at the same time. 1 As a result, we expect that new 
models which have the ti produced via t- or u-channel are subject to similar restrictions 
when fine-tuning for a workable parameter set. 

Similarly, using the definition in Eq. ffl4l) . the ratio Rcp as a function of g' 2 and mw' is 
plotted in Fig. [BJ where (a)-(c) correspond to the decays Bd — > ir + K" , B~ — > ir°K~ and 
Bd — > 7T° K°, respectively. Looking at Fig. E](b), it seems that there is a sizable effect on the 
CPA for B~ — y ii K~ . However, because the SM result for Acp(B~ —y tt°K~) is small, the 



1 Note that the choice of the numerical values used here is for convenience and ease of comparison with 
existing work such as This selection of parameter space would not in any way bias our conclusion. 
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FIG. 4: Contours for the ratio Rb as a function of g' 2 and mw, where the correspondence of each 
plot is (a) B~ -> tt~K°, (b) B d -> vr+A^, (c) B" -> vr°A^ and (d) 5 d -> v^A . 
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FIG. 5: Allowed range (gray) for g' 2 and mw, where the world average BR for B — > tt A with 
2cj errors has been used. 

resulting shift in magnitude of CPA due to the presence of new physics is only around 1%. 
Furthermore, if we combine the result in Fig. [6]^a) with the allowed range depicted in Fig. |5j 
we see that the modification from the W effects will be about 10% of A^p{Bd — > tt + K~) 
only. Therefore, the change in the absolute value of CPA for — > n + K~ is at about the 
1% level. Finally, we note from Fig. |6]^c) that the CPA for Bd — > 7r°K° is more or less 
insensitive to the W effects. 



9 



loo 500 600 700 800 JoO 500 600 700 800 JoO 500 600 700 800 

m W i [GeV] raw 1 [GeV] raw 1 [GeV] 



FIG. 6: Ratio of CPA in the W exchange to that in the SM for (a) B d -> tt+K", (b) B~ -> 
and (c) B d -> ir°K°. 

In summary, the fascinating discovery of an unexpectedly large FBA in the top quark 
pair production at the Tevatron strongly suggests that some new physics may be at play. 
While the exact nature of this is not clear, we observe that any new interactions introduced 
to solve the anomaly which produce tt via the t- or w-channel are naturally correlated to 
rare B decays. In this work, we illustrate this connection by studying a model with a new 
charged gauge boson W added to the SM. Using the current precision measurements on 
B — > tcK decays, we estimate the constraints on the parameters of the theory. In particular, 
we find that strong conditions are imposed on the new gauge coupling and the mass of the 
associated gauge boson. From the data for B~ — > ir~K , which is the most stringent among 
the set of B — > nK decays considered, our study shows that the size of the coupling g' 2 is 
strongly constrained from above for the relevant W mass range of several hundreds GeV. 
On the other hand, these W-mediated effects are not expected to play a significant role on 
the direct CPA of B decays. 
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Appendix A: formulations of / BM (0) and 9o //lA/2 ( m D in PQCD approach 

For calculating the transition and time-like form factors in PQCD approach, the necessary 
distribution amplitude of B meson and outgoing light meson is defined by 

(0Mz)b a (0)\B(p B )) = J d 4 A;e^|(^ + m B ) 75 B (fc) - Mk) J , 



(M(p)\q 2 p(z)q la (0)\0) = --== J dxe*"* { 7 5 &>m{x) + >y 5 m M & M (x) 

+m° A MA + A--mii} a p ■ (Al) 

with p B = 1, 0±),p= (p + , 0, 0_l), n+ = (1, 0, ± ), n_ = (0, 1, 0j_). Since the effects of 

<Pb are small 50|, in the calculations we only focus on the contributions from By hard 
gluon exchange, the B — » M transition form factor is expressed by 

pi POO 

/ SM (0) = %nC F m B / dx x dx 2 \ b 1 db 1 b 2 db 2 <f> B (^lA) 
Jo Jo 

{[(1 + x 2 ) $ M (x 2 ) + r M (l - 2x 2 ) [§ P M (x 2 ) + $^ (x 2 ))} E e 
h e (x u x 2 , b u b 2 ) + 2r M $ p M (x 2 ) E e (t®) h e (x 2 , x u b 2 , 6i)} (A2) 



x 

X 



with Cf = 4/3 and tm = rn^/mB, and the — > M\M 2 time- like form factor is written as 

M\M% t 2n {m b + m u ){m s - m d ) 2 f 1 f°° 

% \™b) = 2 2 SirC F m B / dx 2 dx 3 / b 2 db 2 b 3 db 3 



x { [r M2 x 2 <S> Ml (x 2 ) ($m 2 (1 - x 3 ) + $^ 2 (1 - x 3 )) + 2r Ml $P /i (x 2 )$ Af2 (l - x 3 
xE a (tl) h a {x 2 ,x 3 ,b 2 ,b 3 ) + [x 3 r Ml (® p Ml {x 2 ) - ^° Ml {x 2 )) $m 2 (1 - £3) 
+2r M2 $A/ 1 (x 2 )$P /2 (l - x 3 )] E a (t 2 a ) h a (x 3 , x 2 , 6 3 , & 2 )} • (A3) 
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The hard functions /i e , a are given by 

h e (x 1 ,x 2 ,b 1 ,b 2 ) = S t (x 2 )K (y/x 1 X2m. B b 1 ) 

x[0(h - b 2 )K (^/x^m B b 1 )I (^/x^m B b 2 ) 



+9(b 2 - b 1 )K (^ 2 mBb2)Io{V%^m B bi)) , (A4) 

2 



X 



6(b 2 - b 3 )H^ (yj x 3 m 2 B b 2 )J (^ x 3 m 2 B b 3 

+ 6{b 3 ~ b 2 )H£\^x 3 mlb 3 )J (^rf B b 2 ) 
and the evolution factor E ea are defined by 



(A5) 



E e (t) = a s (t) S B (t) S P (t) , 

E a (t) = a s (t)S Pl (t)S P2 (t), (A6) 
where Suit) denotes the Sudakov factor of M-meson and S t is the threshold resummation 



effect. Their explicit expressions could be found in Ref. 



45] and the references therein. The 



hard scales for emission and annihilation are chosen to be 



t\ = max(y / x 2 m|, l/b u l/b 2 ) , 
t\ = max(Jx 1 m%, 1/h, l/b 2 ) , 



t\ = max(y / x 3 m|, l/b 2 , 1/63) , 

t 2 a = max(Jx 2 m 2 B ,l/b 2 ,l/b 3 ). (A7) 



[1] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 100, 142002 (2008) |arXiv:0712.085Tl 
[hep-ex]]. 

[2] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 101, 202001 (2008) |arXiv:0806.2472l 
[hep-ex]]. 

[3] T. Aaltonen et al. [CDF Collaboration] , |ar2Gvill0L0034J [hep-ex] . 

[4] J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 (1999) [arXiv:hep-ph /9905386] . 
[5] O. Antunano, J. H. Kuhn and G. Rodrigo, Phys. Rev. D 77, 014003 (2008) [arXiv: 0709 .16521 
[hep-ph]]. 

12 



[6] A. Djouadi, G. Moreau, F. Richard and R. K. Singh. la7Xiv:0906.0604l [hep-ph]. 

[7] R Ferrario and G. Rodrigo, Phys. Rev. D 80, 051701 (2009) |arXiv: 0906 .554111 [hep-ph]]. 

[8] S. Jung, H. Murayama, A. Pierce and J. D. Wells, Phys. Rev. D 81, 015004 (2010) 

|arXiv:0907.4TT2l [hep-ph]]. 
[9] K. Cheung, W. Y. Keung and T. C. Yuan, Phys. Lett. B 682, 287 (2009) |arXiv:0908.2589l 

[hep-ph]]. 

[10] P. H. Frampton, J. Shu and K. Wang, Phys. Lett. B 683, 294 (2010) |arXiv:091 1.29551 [hep- 
ph]]. 

[11] J. Shu, T. M. P. Tait and K. Wang, Phys. Rev. D 81, 034012 (2010) |arXiv:0911.3237l [hep-ph]]. 
[12] A. Arhrib, R. Benbrik and C. H. Chen, Phys. Rev. D 82, 034034 (2010) [ arXiv:0911.4875l 
[hep-ph]]. 

[13] P. Ferrario and G. Rodrigo, JHEP 1002, 051 (2010) |arXiv:0912.0687l [hep-ph]]. 

[14] I. Dorsner, S. Fajfer, J. F. Kamenik and N. Kosnik, Phys. Rev. D 81, 055009 (2010) 

|arXiv:0912.0972l [hep-ph]]. 
[15] D. W. Jung, P. Ko, J. S. Lee and S. h. Nam, Phys. Lett. B 691, 238 (2010) [arXiv:0912.1105l 

[hep-ph]]. 

[16] J. Cao, Z. Heng, L. Wu and J. M. Yang, Phys. Rev. D 81, 014016 (2010) [arXiv:0912. 14471 
[hep-ph]]. 

[17] Q. H. Cao, D. McKeen, J. L. Rosner, G. Shaughnessy and C. E. M. Wagner, Phys. Rev. D 
81, 114004 (2010) |arXiv:1003.346T1 [hep-ph]]. 

[18] K. Kumar, W. Shepherd, T. M. P. Tait and R. Vega-Morales, JHEP 1008, 052 (2010) 
|arXiv: 1004.48951 [hep-ph]]. 

[19] M. V. Martynov and A. D. Smirnov, iarXiv: 1006.42461 [hep-ph]. 

[20] R. S. Chivukula, E. H. Simmons and C. P. Yuan. IarXiv: 1007.02601 [hep-ph]. 

[21] M. Bauer, F. Goertz, U. Haisch, T. Pfoh and S. Westhoff. la?Xiv:1008.0742l [hep-ph], 

[22] C. H. Chen, G. Cvetic and C. S. Kim, Phys. Lett. B 694, 393 (2011) |arXiv: 1009.41651 [hep- 
ph]]. 

[23] D. W. Jung, P. Ko and J. S. Lee, larXiv:1011.5976l [hep-ph]. 

[24] E. Alvarez, L. Da Rold and A. Szvnkman. larXiv:1011.6557l [hep-ph], 

[25] D. Choudhury, R. M. Godbole, S. D. Rindani and P. Saha. IarXiv: 1012.47501 [hep-ph], 

[26] K. Cheung and T. C. Yuan. larXiv:1101.1445l [hep-ph]. 

13 



[27] B. Xiao, Y. K. Wang, Z. Q. Zhou and S. h. Zhu, Phys. Rev. D 83, 057503 (2011) 

|arXiv:1101.2507l [hep-ph]]. 
[28] C. Delaunay, O. Gedalia, S. J. Lee, G. Perez and E. Ponton. larXiv:1101.2902l [hep-ph]. 
[29] E. L. Berger, Q. H. Cao, C. R. Chen, C. S. Li and H. Zhang, arXiv:1101.5625l [hep-ph]. 
[30] V. Barger, W. Y. Keung and C. T. Yu, Phys. Rev. D 81, 113009 (2010) j arXiv:1002.1048l 

[hep-ph]]. 

[31] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[32] V. Barger, W. Y. Keung and C. T. Yu, Phys. Lett. B 698, 243 (2011) [ arXiv:1102.0279l 
[hep-ph]]. 

[33] B. Bhattacherjee, S. S. Biswal and D. Ghosh. larXiv:1102.0545l [hep-ph]. 

[34] K. Blum et al. larXiv:1102.3133l [hep-ph]. 

[35] K. M. Patel and P. Sharma. larXiv: 1102.47361 [hep-ph], 

[36] G. Isidori and J. F. Kamenik. larXiv: 1103.00161 [hep-ph]. 

[37] R. Foot, arXiv: 1103. 19401 [hep-ph]. 

[38] C. Delaunay, O. Gedalia, Y. Hochberg, G. Perez and Y. Sorea. larXiv":1103.2297l [hep-ph], 
[39] J. A. Aguilar-Saavedra and M. Perez- Victoria, arXiv:1103.2765 [hep-ph]. 
[40] M. I. Gresham, I. W. Kim and K. M. Zurek. IarXiv:1103.3501l [hep-ph], 
[41] J. Shu, K. Wang and G. Zhu. larXiv:1104.0083l [hep-ph]. 

[42] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 106, 171801 (2011) | arXiv: 1104.06991 
[hep-ex]]. 

[43] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996) 
[arXiv:hep-ph/9512380] . 

[44] H.N. Li and G. Sterman, Nucl. Phys. B381, 129 (1992); G. Sterman, Phys. Lett. B179, 281 
(1986); Nucl. Phys. B281, 310 (1987); S. Catani and L. Trentadue, Nucl. Phys. B327, 323 
(1989); Nucl. Phys. B353, 183 (1991); H.N. Li, Phys. Rev. D64, 014019 (2001); H.N. Li, Phys. 
Rev. D66, 094010 (2002). 

[45] C. H. Chen, Y. Y. Keum and H. n. Li, Phys. Rev. D 66, 054013 (2002) |arXiv:hep-ph/0204166] . 

[46] C. H. Chen, C. S. Kim and Y. W. Yoon, Phys. Lett. B 671, 250 (2009) [ arXiv:0801.0895l 
[hep-ph]]. 

[47] H. n. Li, S. Mishima and A. I. Sanda, Phys. Rev. D 72, 114005 (2005) [arXiv:hep-ph/0508041] . 
[48] The Heavy Flavor Averaging Group et al., larXiv: 1010. 1 589 [hep-ex]. 



11 



[49] H. C. Yen and I. F. Shih, Chin. J. Phys.13, 205 (1975). 

[50] C. D. Lu and M. Z. Yang, Eur. Phys. J. C 28, 515 (2003) [arXiv:hep-ph/0212373] . 



15 



